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ABSTRACT 

We use the most up-to-date cosmological evolution models of dusty star-forming galaxies and radio sources to compute the extra- 
galactic number counts and the cosmic background from 408 MHz to 12 THz (or 24 /jm). The evolution model of star-forming 
galaxies reproduces the recent constraints obtained by Spitzer, Herschel, and ground-based submm and mm experiments: number 
counts, redshift distribution of galaxies, cosmic background intensity and anisotropies. The template spectral energy distributions 
used in this model are extrapolated to the radio domain adding three emission components: synchrotron, free-free, and spinning dust. 
To fix the synchrotron intensity, we use the well-known IR/radio flux ratio, qio, and a constant spectral index /Ss = -3, consistent 
with measurements made in local galaxies taking account the spinning dust emission. For a constant q 10 , our model added to the 
AGN contribution provides a good fit to the extragalactic number counts from 24 fjm to 408 MHz, and to the cosmic background 
intensity in the far- and mid-IR. The spinning dust emission accounts for up to 20% of the cosmic microwave background produced 
by star-forming galaxies, but for only less than 10% of the total background when AGN are included. The star-forming galaxies 
account for 77.5% of the number counts at 1.4 GHz for a flux of 100 yujy. However, the model falls short of reproducing the cosmic 
radio background measured with the ARCADE2 balloon-borne experiment. Considering the case when qjo decreases strongly with 
redshift, this still does not explain the ARCADE2 measurements. It also yields to an overestimate of the low-flux number counts in 
the radio. As a result, we rule out a steep variation of qio with the redshift at least for z < 3.5. Then, adding a population of faint 
star-forming galaxies at high redshift (L IR < 10" L a and 4 < z < 6), which would be able to reproduce the ARCADE2 measurements, 
leads to predictions of the cosmic IR background much higher than what is currently observed, ruling out this as the explanation for 
the ARCADE2 results. Considering our findings and previous studies of the diffuse extragalactic radio emission, we conclude that if 
the radio emission measured by ARCADE2 is astrophysical in origin, it has to originate in the Galaxy or to originate in a new kind of 
radio sources (with no mid- to far-IR counterparts) or emission mechanism still to be discovered. 

Key words. Galaxies: evolution - Cosmic background radiation - Infrared: galaxies - Radio continuum: galaxies 



1. Introduction 



Recent measurements of the CRB made by iFixsen et al.l 



The extragalactic background light or cosmic background is 
dominated by the Cosmic Microwave Background (CMB), a 
relic of the evolution of the early Universe. The remainder of 
the background is directly linked to the formation and evolu- 
tion of all the galaxies in the Universe apart from our Galaxy. 
Even if less energetic, this background can teach us a lot about 
galaxy evolution throughout time and can be separated into five 
components: the X- and y-ray backgrounds, the cosmic op- 
tical background, the Cosmic Infrared Background (CIB), and 
the Cosmic Radio Background (CRB). The CIB is produced 
by the emission of dust illuminated by the stars contained in 
infrared (IR) or dusty star-forming galaxies. The evolution of 
these galaxies is well constrained by direct measurements of 
the CIB, galaxy number counts from the mid- to the far-IR, 
and redshift distribution of galaxies. Star-forming galaxies are 
also expected to emit at radio frequencies because the stars they 
contain can ionise their gas content and thus produce free-free 
emission. Eventually, these stars will evolve into supernovae 
leading to synchrotron emission. Anomalous microwave emis- 
sion (AME), assumed to be produ ced by spinning dust grains 
dPlanck Collaboration et al.|l201 ldllah . may also arise from dusty 
star-forming galaxies. The exact intensity and evolution of radio 
emission from star-forming galaxies is still a matter of debate. 



(I201 lh with the ARCADE2 experiment have renewed the in- 
terest in the contributi on of star-forming g alaxies at microwave 
wavelengths. Indeed, IFixsen et al.l d201 lh measured a back- 
ground between 3 and 10 GHz that is at least 5<x above the 
Cosmic Microwave Background (CMB) te mperature es ti mated 
by COBE/FIRAS. Investigations made by ISingal etal] d2010h 
showed that standard extragalactic radio sources cannot account 
for this radio excess and they concluded that it has to come 
from star-forming galaxies. However, up till now, no extragalac- 
tic source evolution model has been able to explain this excess 
dMassardi et al.ll201Ct ISingal et alj|20lol: IVernstrom et al.l l201 U 
iPonente et al.l [Draper et al.ll201 ll) . Our aim is to provide 
a comprehensive evolution model of star-forming galaxies able 
to reproduce all the observations, cosmic background intensity 
and extragalactic number counts, from the mid-IR to the radio. 
To do so, we use the most up -to-date galaxy evolut ion models 
for both sta r-forming galaxies dBethermin et al.l201 1 ) and radio- 



forboth star-forming galaxies (Bethermin et al. 201 1 
loud AGN d~de Zotti et al.l2005tlMassardi et al.l201oT 



Tucci et all 



1201 lh . Then, we test several hypothesis to try to explain the ex- 
cess measured by the ARCADE2 team: extragalactic spinning 
dust emission? Variation of the well-known far-IR/radio flux ra- 
tio with redshift? Additional population of faint star-forming 
galaxies at high redshift? 

The paper is organised as follows. Section [2] describes the 
evolution model of dusty star-forming galaxies we use, as well as 
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the emission components added to explain radio emission (syn- 
chrotron, free-free, spinning dust). The models for the radio-loud 
AGN are described in Sect. [3] Our results are then presented in 
Sect. [4] Finally we give our conclusions in Sect. [5] 



2. Evolution model for dusty star-forming galaxies 

2.1. Luminosity function evolution and IR spectral energy 
distributions 



We used the recent model of Bet hermin et all d201 ll) to compute 
the number counts of dusty star-forming galaxies at the radio 
wavelengths. Two ingredients come into play: the spectral en- 
ergy distribution (SED) of galaxies and the luminosity function 
(LF) evolution. 

The SEDs, from 4 to 1 mm, are from the template library 
of lLagache et al.l (120041) . which consists of two galaxy popu- 
lations. First, a star-forming galaxy population with SEDs that 
vary with IR bolometric luminosities. Second, a normal-galaxy 
population with a fixed template SED that is colder than the star- 
forming galaxy templates. The normal and star-forming galaxies 
are dominant at low- and high-luminosity, respectively. The frac- 
tion of each galaxy population as a function of the bolometric 
luminosity is given by a smooth function 



tstarburst I + tll[log m (L lR / L pop ) / CT pop \ 



(1) 



where th is the hyperbolic tangent function, L pop the luminosity 
at which the number of normal and star-forming galaxies are 
equal, and cr pop characterises the width of the transition between 
the two populations. At Lir 

50 %. 

IBethermin et al.l (1201 ll) assume that the luminosity function 
(LF) is a classical double exponential function: 



L pop , the starbursts fraction is 



(2) 



where ^(Lir) is the number of sources per logarithm of lumi- 
nosity, and per comoving volume unit for an infrared bolomet- 
ric luminosity Lir. <t>* is the normalisation constant character- 
ising the density of sources, L* is the characteristic luminos- 
ity at the break, and I - a and 1 - a - l/<x 2 /Zn 2 (10) are the 
slopes of the asymptotic power-law behaviour at low- and high- 
lumino sity respectively. To rep roduce the observations in the in- 
frared, IBethermin et al.l d201 ll) allow the parameters of the LF 
to vary with redshift. In order to limit the number of free pa- 
rameters and to avoid strong degeneracies, a and cr are fitted 
as constants with z. A continuous LF redshift-evolution in lu- 
minosity and density is assumed following broken power-laws: 
L* oc (1 + zY L and <I>* oc (1 +zY' c , where ri and are parameters 
driving the evolution in luminosity and density, respectively. It is 
imposs ible to reproduce the ev olution of the LF with constant ri 
and r*. IBethermin et all ( 1201 lb consequently allow their value to 
change at two specific redshifts: three different values of ri and 
are used at low redshift (0 < z < Zbreak), intermediate redshift 
(Zbreak < z < 2), and high redshift (z > 2). The position of the 
first redshift break is a free parameter and converges to the same 
final value (Zbreak ~ 0.9) for initial values < z < 2. The second 
break does not appear naturally in the fitting procedure and is 
fixed at z — ^3- 



The model has 13 free parameters that were determined by 
fitting the model to published measurements of galaxy number 
counts and monochromatic LF measured at given redshifts: 

- Number counts: Spitzer counts at 24, 70 and 160 pm, 
Herschel counts at 250, 350 and 500 pm, and AzTEC counts 
at 1.1 mm. 

- Monochromatic LFs: IRAS local LF at 60 pm, Spitzer LF at 
24 pm at z=0, at 15 pm at z=0.6, at 12 pm at z=l , and at 
8 pm at z=2. 

- FIRAS CIB spectrum between 200 pm and 2 mm. 

The best-fit parameters, as well as their uncertainties and de- 
generacies, were obtained using a Monte Carlo Markov chain 
(MCMC) Metropolis-Hastings algorithm. The model is adjusted 
to reproduce deep counts and monochromatic LFs at key wave- 
lengths. It also reproduces rece nt observations that constrain 
very well the models, such as the Ijauzac et al.l (1201 ll) measured 
redshift distribution of the CIB. 

We used the so-called mean model, which i s obtained using 
the m ean value of the parameters as given in Beth ermin et al.l 
(120111) with the lensing contribution that is available on the 
http : //www. ias .u-psud. f r/irgalaxies / web p age. 

The SED templa tes of lLagache et al.l (120041) used in 
IBethermin et al.l (1201 ll) stop at A — 1 mm. To predict the num- 
ber counts and cosmic background at radio wavelengths of dusty 
star-forming galaxies, we need to extrapolate the SED templates 
into the radio domain. We consider three emission components: 
synchrotron (Sect. 12.21 ). free-free (Sect. 12.31 ). and spinning dust 
(Sect. Ell. 

2.2. Synchrotron 

Synchrotron emission is produced by relativistic electrons inter- 
acting with a magnetic field. At first order, its emissivity can 
be expressed as a power-l aw: es (v) °c v 85 . Using WMAP inten- 
sity and polarisation data. lMiville-Deschenes et aD (120081) mea- 
sured a spectral index [}$ = -3.0 + 0.06 in our Galaxy between 
408 MH z and 23 GHz, in ag reem ent with previous m easure- 



ments bv lGiardino et alJ (120021) andlPlatania et aD(l2003h . In the 
nearby spiral galaxy M33. lTabatabaei et al.l d2007l) showed that, 



1 Numerical values and details about the fitting procedure are given 
in lBefhermin et all |20TH) . 



even if the synchrotron spectrum hardens from the outer parts 
of the galaxy to the star-forming regions, a constant emissivity 
spectral index, /3s = -2.953 ± 0.285 , is appropriate fo r global 
studies. Similar values were found by iPeel et al.l d201 ll) in three 
nearby star-forming galaxies (M82, NGC253, and NGC4945). 
Consequently, and for the sake of simplic ity, we adopt a single 
power-law for the two galaxy templates of lLagache et al.l (120031) 
with y8 s = -3.0. 

The radio continuum and far-IR emi ssions are observe d to be 
tightly correlated in external galaxies dHelou et alj|1985l ). This 
correlation is remarkable in the sense that it holds over five 
order s of magnitude of radio continuum luminosity dYun et al.1 
l200lh . an d both at the scale of an entire galaxy and at subgalac- 
tic scales dTabatabaei et alJl2007t iHughes et al]l2006l). The cor- 
relati on is also observed in the Milky Way dBoulanger & Perault 
| 1988[) a t large scale and down t o a few parse cs scale dZhang et al. 
bold) . iBoulanger & Peraultl dl988l) and IZhangetalJ (l2010h 
found that the far-IR/radio flux ratio, usually called q, has sim- 
ilar values in the Milky Way and in external galaxies. Because 
synchrotron is the dominant emission component at 1.4 GHz, 
we set its intensity using the far-IR/radio flux ratio, qw = 
logiSjQ/tm/S i.4GHzX of IR emission at 70 pm and radio emission 
at 1.4 GHz. 
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The possible evolution of ^70 with reds hift is still de- 
bated, some authors finding a constant relation (llbar et al 1120081 
ISargent et al.l |2010|). where others find a slow decrease with 
redshift dSevmour et al.ll2009t iBeswick et al.ll2Q08i: llvison et all 
I2010bl) . As a result, we will consider two cases in this paper. 
The first one assumes that the far-IR/radi o flux ratio is constant 
over redshift and equal to qm ~ 2.15 (lAppleton et al.l 120041: 
iFraver et al.l 120061: ISevmour et al.ll2009l) . In the following, we 
will refer to it as case A. The resulting spectra for the nor- 
mal spiral galaxy and starbust galaxy templates are displayed in 
Fig- HI The second one assumes that the ratio evolves with red- 
shift as qjoiz) ~ 2.15 — 0.72 log(l + z). This relation comes from 
ISevmour et al.l (l2009h and is the steepest that has been reported 
so far. We choose it to maximize the effect of varying qjo on ex- 
tragalactic number counts. We will refer to it as case B. Case B 
will lead to a brighter synchrotron emission for high redshift and 
should consequently have a stronger impact on the prediction of 
the CRB than case A. 



2.3. Free-free 

The free-free component is particularly difficult to e stimate as 
it is ne ver dominant at any wave length. According to lPeel et al.l 
d2011l) and lMurphv et al.l ( 1201 Ol) . we model the free-free spec- 
trum as a power-law with a spectral index Bef = -2.1. 
This value is similar to what is measured in ou r Galaxy 
dMiville-Deschenes et"aill2008l:lDickinson et~ai1l2003l) . In previ- 
ous studies, before anomalous microwave emission (or spinning 
dust) was observed in galaxies, the free-free intensity was of- 
ten norm alised to be equal to that of the syn chrotron at about 
30 GHz dCondonlll992T iPonente et ail 1201 ll) . This means that 
spinning dust emission was included in the free-free and/or syn- 
chrotron emission. However, the spectral shape of the spinning 
dust emission is a "bump" centred around 20-30 GHz while the 
free-free and synchrotron emissions span over a much broader 
range of frequencies. Thus, ignoring the spinning dust contribu- 
tion can lead to a flatter spectral index for the free-free and/or 
synchrotron emissions that implies an overestimate of their in- 
tensity. Consequently, we use recent results taking into account 
spinning dust emission to normalise the free-free emission at 
10 GHz. For t he thr ee nearby star-forming galaxies observed 
by iPeel et alJ (1201 lh . the free-free emission is about half of 
that of the sync hrotron at 10 GHz. S imilar values are found 
in NGC6946 by iMurphv et al] d2010l) : the free-free intensity 
reaches about (52 + 12)% of the synchroton intensity at 10 GHz 
th roughout the galaxy. Th ese values also match the observations 
of iTabatabaei et all d2007l) forM33. Consequently, we normalise 
the intensity of the free-free emission to be equal to 50% of the 
synchrotron at 10 GHz. The resulting spectra for the normal spi- 
ral galaxy and star-forming galaxy templates are displayed in 
Fig.[TJfor case A. 



2.4. Spinning dust emission 



microwave emission, peaking between 10 and 
is observed everywhere in our Galaxy, from 



Anomalous 
60 GHz, 

the Galactic plane dPlanck Collaboration et al. 1201 lal) to in- 
termediate G alacti c latitudes JMivS-Deschenes et al' l r2003 
iDavies et al. 1 120061: lEagachel 120031) . AME is detected in var- 
ious interstellar clouds from the most diffuse to dense 
molecular clouds dPlanck Collaboration et al.l201 ldt Vidal et al. 
l201lt ICasassus et all 120081 120061: IWatson et al.l 120051) . Up to 
now, a few extragalactic detections have been reported in 



nearby star-forming gal axies dPeel et alJ 1201 it IMurphv et al.l 
12010b IScaifeet al.l 1201 Oh and in the Small Magel lanic Cloud 
dPlanck Collaboration et ail 1201 Id iBot et all l2Q10t) . It is usu- 
ally assumed that AME is due to the emission of rapidly ro- 
tating interstellar P AHs, the so-called spinnin g dust emission, 
as first suggested by iDraine & Lazarianl (119981) . The rotation of 
the PAHs is driven by interactions with the interstellar gas par- 
ticles and absorptions of UV/visible stellar photons. This means 
that any environment containing PAHs should emit in the mi- 
crowave. The spectra of the dusty star-forming galaxies respon- 
sible for the CIB exhibit bands in the mid-IR (~ 3 to 20 fim) 
characteristic of the interstellar PAHs. Consequently, we add a 
spinning dust component to our galaxy templates to estimate its 
impact on the source counts and cosmic background intensity in 
the microwave and radio wavelengths. 

The intensity and the peak frequency of the spinning dust 
spectrum depend on the strength of the interstellar radiation field 
heating the grains and on the local gas dens ity and temperature 
dDraine & Lazarianll998blYsard et al.l201 ll) . Thus, two parame- 
ters have to be set to determine the spinning dust component for 
star-forming galaxies: a mixture of environments and the radia- 
tion field intensity. 

For the Milky Way, the comp onent separation performed 
by iMiville-Deschenes et al.l (120081) using WMAP data allowed 
I Ysard et alJ (l2010h to measure a ratio between the AME in the 
WMAP 23 GHz band and the PAH emission in the IRAS 12 fim 
band of L 2 T,GHz/Ln f im = (5.36 x 10~ 3 ) ± (4 x 10~ 5 ) for the en- 
tire Galaxy. This ratio is valid for the radiation field illuminating 
our Galaxy. The intensity of the radiation field can be expressed 
in units of Go-factor, equal to 1 in the case of the interstel- 
lar standard radiatio n field (ISRF) illuminati ng t he solar neigh- 
bourgh ood for which iBoulanger et alJ d 1996b and lLagache et alJ 
d 1998b measured an average dust temperature of ~ 17.5 K. 
Using the FIRAS spectrum and the DIRBE 100 fim data aver- 
aged over the entire Galaxy, we determine the mean dust tem- 
perature of the Milky Way by fitting a modified blackbody with 
a constant spectral index B = 2 for 100 ^ A < 250 yum. 
We find a dust temperature equal to Tsed ~ 18 K leading to 
G mw ^ (T SED /n.5 K) 4+ ^ = 1.2 (lYsard et alJl20Toh . Then, we 
proceed to the same fitting procedure to get estimates of the in- 
tensities of the mean radiation fields required to explain the far- 
IR spectra of our galaxy templates. For normal spiral galaxies, 
we get Go ~ 2 {Tsed = 19.5 K), and for star-forming galaxies 
with bolometric luminosities of L« = 10 10 , 10 11 , and 10 12 L , 
we get Go ~ 13, 22, and 45, respectively (Tsed ~ 27, 29, and 
33 K, respectively). Note that the exact B value and wavelength 
range used to compute the temperature are not critical as long as 
the same parameters are used for both the Galactic spectrum (the 
reference) and the extragalactic SED templates. 

Then, the spectral shape of the spinning dust sp ectrum is de- 
termin ed with the spinning dust model described in lSilsbee et al.l 
(1201 ll) . We assume that the ISM of galaxies can be approximated 
with a mixture of environments dDopita & Sutherlan dl2003t lCoxl 
120051) composed by: 20% of Warm Ionised Medium (WIM, with 
a proton density of = 0.1 cm -3 ), 20% of Warm Neutral 
Medium (WNM, n H = 0.4 cm 4 ), 30% of Cold Neutral Medium 
(CNM, n H = 30 cm -3 ), and 30% of Molecular Cloud (MC, 
tin = 300 crrT 3 ), and illuminated by the ISRF scaled with the 
Go values determined previously. Then, we integrate the galaxy 
template spectra in the IRAS 12 fim filter in order to normalise 
them to the L23GHz/^i2^m ratio measured bv lYsard et al.l (1201 Ol) 
in the Milky Way. To take into account the lower Gq* 1 * = 1.2 
value measured in the Milky Way, we multiply the spinning dust 
spectrum by a factor equal to the ratio of the integrated inten- 
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Fig. 1. SED template for a galaxy with a bolometric luminosity 
of 10" L Q (black line), contribution of IR dust emission (green 
line), free-free emission (blue line), synchrotron emission cal- 
culated for case A (cyan line), and spinning dust emission (red 
line), a) Normal spiral galaxy template, b) Star-forming galaxy 
template. 



Table 1. References for the extragalactic number counts pre- 
sented in Figs. |2 3] and [5] 

Frequency References 



408 MHz 
610 MHz 
1.4 GHz 

2.7 GHz 

4.8 GHz 
5 GHz 
8.4 GHz 
15.7 GHz 
20 GHz 
100 GHz 
143 GHz 
145 GHz 
148 GHz 
214 GHz 
217 GHz 
250 yum 
24fim 



Gervasi et al. (2008a) and references therein 
^ervas^et^ah (2008a) and references therein 
Gervasi et al. (2008a) and references therein 
^ervas^et^aL (2008a) and references therein 
Vernstrom et al. (2011) 
^lervasrirt^ah (2008a) and references therein 
Tuccietal. (2011) 
AMI Consortium et al. (2011) 
Tuccietal. (2011) 
F^lanckCohaboration et_al. 
Planck Co llaboratio n et al 
Vieira et al. (2010) 
Marriage et al. (2011) 
Vieira et al. (2010) 
r^lanckC^llaboration^tjiL 
Bethermin et al. (2012) 
Bethermin et al. (2010) 



sity in the WMAP 23 GHz band of this spinning dust spectrum 
computed for a given Go value and the same but computed for 
qmw -12. This factor is equal to 1 .00 in the case of normal spi- 
ral galaxies, and to 1.15, 1.28, and 1 .45 for star-forming galaxies 
with bolometric luminosities of Ljr = 10 10 , 10 11 , and 10 12 L Q , 
respectively. 

The resulting spectra are displayed in Fig.[T]for case A (con- 
stant qjo). For normal spiral galaxies, the spinning dust emission 
at 23 GHz accounts for 31.7% of the total intensity. For star- 
forming galaxies with L/r = 10 10 , 10 n , and 10 12 L , it accounts 
for 25.3%, 18.5%, and 12.8% of the total intensity, respectively. 
If qjo varies with the redshift (case B), the spinning dust com- 
ponent is then fully dominated by the synchrotron and free-free 
emissions. 



3. Evolution model for the AGN-powered radio 
sources 

The extragalactic radio number counts are dominated by radio- 
loud ac tive galactic nuclei (AGN) for b right flux densities, S > 
1 mJy dlbar et al.ll2008l: ICondonlll992l) . These radio-loud AGN 
are also expected to account for about half of the surface bright- 
ness of the CRB dMassardi et al.ll2010l) . To take into account the 
contribution of radio galaxies, we use two different models, de- 
pending on the frequency. For v = 408 MHz, 6 1 MHz, 1 .4 GHz, 
2.7 GHz, 5 GHz, 1 5.7 GHz, 30 GHz, 44 GHz, and 70 GHz, 
we use the model of iMassardi et al.l d2010l) . For v = 8.44 GHz, 
20 GHz 100 GHz, 1 43 GH z, and 214 GHz, we use the model 
C2Ex of lTucci et al.l j20Tlh . 



These two models provide good fits to the extragalactic num- 
ber counts in the radio but the model of lTucci et al.l (1201 ll) gives 
a better repartition of the coun ts between the flat- a nd steep- 
spectrum galaxies. The model of Mass ardi et alj ( 120101) assumes 
two flat-spectrum populations with different luminosity func- 
tions and one steep-spectrum population. For each population, 
they adopt a power-la w spectrum S v cc v~" with a = 0. 1 or 
a = 0.8. The model of lTucci et al.l (1201 ll) allows for a gaussian 
distribution for the spectral indices. They also introduce an addi- 
tional population with an inverted spectrum and define frequency 
breaks for the slopes of the flat-spectrum populations. 




0.00010 
0.00001 



10 100 1000 

Frequency (GHz) 



Fig. 3. Cosmic background spectrum predicted by our model 
for dusty star-forming and normal spiral galaxies without (dot- 
ted line) and with (dashed line) spinning dust emission. The 
solid line shows the s pectrum when the contribution of radio- 
loud AGN is inclu ded dde Zotti et alj|2005l: IMassardi et al.120101 
iTucci et aill201 lb . The blue triangles sh ow direct measurem ents 
of the CRB obtained by ARCADE2 (iFixsen et all l201ll) and 
older surveys reprocessed using the model developed by the 
ARCADE2 team to measure the CRB (408 MHz and 1.4 GHz). 
The green triangles show direct measurements of the CIB at 
250, 350, and 500 //m using FIRAS data processed as in 
lLagache et al.1 d 1999b . The r ed triangle shows the C IB measured 
in the MIPS 24 jum band bv lBethermin et al.l (l2010t) . 



4. Results and discussion 

4.1. Constant q 7 o and spinning dust emission (case A) 

The predictions for the differential extragalactic number counts 
in case A (constant qjo, Sect. I2. 2b are presented in Fig. [2] The 
contributi ons of both dusty star-forming galaxies and AGN are 
includ ed dde Zotti et al.ll2005b IMassardi et al.ll2010l: ITucci et al.l 
1201 ll) . For the dusty star-forming galaxies, Fig. [2] also presents 
the uncertaint ies on the count s calcu lated from the MCMC fit 
of the LF by iBethermin et al.l (1201 ll) in the case of the SEDs 

described in Sects. 12771 12. 21 and l2.4J 

In agreement with IMassardi et al.l (1201 Oh . we find that the 
low-flux bump of the counts in the radio (610 MHz and 
1 .4 GHz) can be accounted for by standard star-forming galax- 
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Table 2. Surface brightness of the cosmic background for case A 
(first column), contribution of AGN (second column), and con- 
tribution of the IR galaxies (third column). 



Frequency 


Total 


AGN 


IR galaxies 


(GHz) 


(nW/nr/sr) 


(nW/nr/sr) 


(nW/m 2 /sr) 


0.408 


6.83X10- 5 


4.16X10- 5 


2.67X10- 5 


0.61 


7.32X10- 5 


4.60x1 0- 5 


2.73X10- 5 


1.4 


9.19xl0" 5 


6.24xl0~ 5 


2.96xl0- 5 


5 


1.25xl0~ 4 


8.51xl0 _5 


3.99xl0~ 5 


8.4 


1.58xl0~ 4 


1.04X10- 4 
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ies. Moreover, our model can explain all the galaxy number 
counts from the mid-IR (24 //m) to the radio (408 MHz). Our 
model shows that star-forming galaxies account for 84.5% of 
the extragalactic number counts at 1 .4 GHz for a flux of 50 juJy 
and 77.5% for 100 (Ay. 

The CRB predicted by our model can be calculated by 
summing the galaxy number counts. The results are shown in 
Fig. [3] and the corresponding surface brightness from 408 MHz 
to 217 GHz are listed in Tab. 14.11 The contribution of star- 
forming galaxies to the CRB reaches 39.1% at 408 MHz, 37.2% 
at 610 MHz, and 32.2% at 1.4 GHz. These val ues are much 
higher than those found by iPonente et alJ (1201 ll) . who used a 
model based on an empirical relation between the 1 .4 GHz lu- 
minosity and the star formation rate density (SFRD). This dis- 
crepancy could be caused by the fact that these authors used a 
SFRD/radio conversion ratio that is the same for all redshifts. 
This ratio should rather be weighted by the luminosity of the 
galaxies that dominate the averag ed SFRD at a given redshift. 
According to lLagache et al. (2004), the IR output energy is dom- 
inated by galaxies with bolometric luminosities around L« ~ 
8 x 10 10 L Q at z — 0, whereas it is domin ated by galaxi e s with 
L IR ~ 3 x 10 12 L G at z = 2. Furthermore. IPonente et all d2011l) 
do not consider spinning dust emission and assume slopes of 
B = -2.7 for the radio emission, respectively. For local galax- 
ies, the measurements made taking into account the presence of 
spinning dust emission show that it is steeper with B 3. 

As can be seen from Fig.[2j the contribution of spinning dust 
to the galaxy number counts is significant at 8.44, 15.7, 20, and 
30 GHz. However, deeper surveys would be required to validate 
the model in this frequency range as the star-forming galaxies 
dominate the counts only for S ~ 10 juJy, where no data are 
available. Then, spinning dust emission accounts for about 10 to 
20% of the cosmic background produced by star-forming galax- 
ies in the microwave (14.8% at 8.44 GHz, 21.6% at 15.7 GHz, 
19.8% at 20 GHz, and 11.1% at 30 GHz). However, this repre- 
sents only less than 10% of the total background when AGN are 
included (5.1% at 8.44 GHz, 8.5% at 15.7 GHz, 8.3% at 20 GHz, 
and 5% at 30 GHz). 

Finally, Fig.[3]shows that even if our model is able to repro- 
duce all the galaxy number counts and the cosmic background in 
the mid- and far-IR, it c annot explain the e xcess radio emission 
reported by ARCADE2 dFixsen et al]l201 ll) . 




10" ! 10"' 
Flux (Jy) 



100.0 




10"' 



Fig. 4. Differential extragalactic number counts (green triangles) 
for v = 610 MHz (a) and 1.4 GHz (b), in the case of a variable 
far-IR/radio flux ratio qio(z) = 2.15 — 0.72 log(l + z) (see case B 
in Sect.l2~2b. 



4.2. Evolution ofq 10 with redshift (case B) 

The predictions of the differential extragalactic number counts 
for case B (qjo decreasing with redshift, see Sect. \2.2\ are pre- 
sented in Fig.|4]at 610 MHz and 1.4 GHz. As expected, allow- 
ing qjo to vary increases the radio emissivity of star-forming 
galaxies and leads to an increase in the CRB of about 70%, 
which is however still not enough to explain ARCADE2 data: 
ARCADE measurements are still greater than the model pre- 
dictions by a factor of 5 at 408 MHz and of 11.4 at 1.4 GHz. 
Then, case B yields an overestimate of the galaxy counts in the 
radio range (Fig. 0}. This suggests that the far-IR/radio flux ra- 
tio should rather be constant or at least decrease less rapidly 
with redshift as suggested by llvison et"ai1 (1201 Obi) . Acco rding 
to the compilation of data made by Ponente et alJ ( 1201 lb . see 
their Fig. 1, it seems reasonable to consider th at the far- 
IR/radio ratio i s constant at first or de r for z < 3.5 (|Belll | 2003 : 
i Murphvl [20091: llvison et alJ l201QalTbi: iMichalowski et al.l bold 
Isargent et al1l2010tlBourne et al.ll201li) ~ 



We test the case of star-forming galaxies having the proper- 
ties of case A for z < 4 and an enhanced radio emission above. 
The value of qjo(z > 4) needs to be equal to ~ for the model 
predictions to reach the ARCADE2 radio excess. However, this 
again yields to an overestimate of the galaxy counts in the ra- 
dio range, even worse than for case B. As a result, even if the 
far-IR/radio flux ratio decreases only for high redshift, it cannot 
explain ARCADE2 excess without being in contradiction with 
extragalactic number counts. 
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4.3. An additional population of faint galaxies at high 
redshift? 

If we suppose that the radio excess observed by ARCADE2 orig- 
inates in star-forming galaxies, then the only remaining way to 
increase the CRB is to consider an additional population of faint 
IR galaxies at high redshift. It is important to consider faint and 
not bright galaxies not to disrupt the galaxy number counts and 
also because we do not expect a large amount of bright galaxies 
to be undetected even at high redshift. 

We use the galaxy templates o f case A (Sect. \ 2 . 2\ and start 
with the luminosity functions of Be thermin et al. ( 201 ll) . We 
multiply them by a constant for 4 < z < 6 and for L/« < L t h res . 
The lowest L,/„ f , allowing to reproduce the ARCADE2 radio 
excess is L,i ires = 10 11 L . For lowest threshold luminosi- 
ties, it is not possible to fit simultaneously the background sur- 
face brightness at 408 MHz and 5 GHz. Higher threshold lu- 
minosities predict too many galaxy number counts at low flux 
(< 1 mJy) that we ruled out by current observational constraints. 
The model predictions for L,i ires = 10 11 L Q are presented in Fig. 
|5] The increase in the luminosity function required to match the 
ARCADE2 excess overproduces the CIB by a factor of 4.6 at 
24 /mi, 8.5 at 250 yum, 1 1.6 at 350 yum, and 21.5 at 500 /im. We 
conclude that a population of faint galaxies at high redshift able 
to account for the radio excess reported by ARCADE2 is ruled 
out by direct measurements of the CIB. 





5. Conclusions 

The aim of this paper was to check if star-forming g alaxies could 
be res ponsible for the radio excess reported by iFixsen et alJ 
d2011l) in the CRB with ARCADE2 data. To do so, using up- 
to-date detailed star-forming and radio galaxy evolution models, 
we considered three scenarios: (/) contribution of spinning dust 
emission, (ii) variation with redshift of the far-IR/radio flux ra- 
tio, and (Hi) additional population of faint IR galaxies at high 
redshift. 

(0 Our model, including the contribution of synchrotron, 
free-free, and spinning dust emissions, and considering a con- 
stant far-IR/radio flux ratio, accounts for all the extragalac- 
tic number counts from the mid-IR to the radio (24 yum to 
408 MHz). It also accounts for the cosmic background in the 
mid- and far-IR. However, it does not reproduce the ARCADE2 
excess. We n o te tha t our model agrees wi th the findings of 
IZannoni et ail d2008l) . iGervasi et aiTd2008bl) . and iTartari et al] 
(120081) . who could explain the TRIS observations with no need 
for an excess in the radio range (600 MHz, 820 MHz, and 2.5 
GHz). Ou r model al s o acco unts for the CRB intensity as esti- 
mated by IWall et all d 19701) at 320 and 707 MHz. Finally, we 
confi rm that radio-lo ud AGN dominate the counts for S v ~Z I 
mJy dlbar et al.ll2008b . 

(ii) In order to increase the radio emission of the star-forming 
galaxies, we considered a possible decrease in the far-IR/radio 
flux ratio, qjo, with redshift. The steepest relation available in 
the litterature does not yield an acceptable fit to the extragalactic 
number counts at 610 MHz and 1.4 GHz. We conclude that vary- 
ing qjo with redshift cannot explain the ARCADE2 radio excess. 
These two results suggest that the far-IR/radio flux ratio should 
rather be constant, or at lea st decrease less rapidly with redshift 
than what was measured bv lSevmour et al.l d2009h . 

(Hi) Finally, we tested the hypothesis of an additional popula- 
tion of faint star-forming galaxies at high redshift. We show that 
the amount of faint galaxies required to explain the ARCADE2 
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Fig. 5. (a), (b), (c) Differential extragalactic number counts at 1 .4 
GHz, 250 pm, and 24 yum (linestyles and symbols are the same 
as in Fig. 13 predicted by our model when a population of faint 
galaxies is added at high redshift (see Sect. I4.3l for details), (d) 
Corresponding cosmic background. Linestyles are the same as 
in Fig. |3] 
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radio excess leads to a predicted CIB much higher than what is 
actually observed at 250, 350, and 500 jjm. 

Consequently, we exclude the star-forming ga laxies as the 
origin of the ARCADE2 radio excess. Furthermore. lSingal et al.l 
d2010t) investigated the contribution to the CRB of radio-quiet 
quasars, radio supernovae, emission from hot gas in galaxy 
clusters, and diffuse emission from regions far from galax- 
ies. They found that they cannot account for more than about 
10% of the excess . Considering our results, and the results of 
ISingal et al.ld20Toh. we conclude that if the radio excess reported 
by Fixsen et al.l ( 201 ll) is astrophysical, it has to have a Galactic 
origin or to originate in a new kind of radio sources (with no 
mid- to far-IR counterparts) or emission mechanism still to be 
discovered. 
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Fig. 2. Differential extragalactic number counts (green triangles) for v = 408 MHz (a), 610 MHz (b), 1.4 GHz (c), 2.7 GHz (d), 
4.8-5 GHz (e), 8.4 GHz (f), 15 GHz (g), 20 GHz (h), 30 GHz (i), 44 GHz (j), 100 GHz (k), 143-145-148 GHz (1), 214-217 GHz 
(m), 250 /im (n), and 24 /im (o). The counts are a compilation from several surveys (see Tab. 14. II for references). The black solid 
and dashed lines are the predictions of our evolution model with and without spinning dust emission, resp ectively, in the case of a 
constant far-IR/radio flux ratio (see case A of Sect. I2.21 > and using the LF-parame ters of the mean model o f lBethermin et al.l ( f201 lb . 
The grey areas show the 1-cr error calculated from the MCMC fit of the L F by [ Bethermin et al. ( 201 1[). The dotted lines are the 



predi ctions of evolution models for radio-loud AGN as described in Sect. [3l (lde Zotti et alJl2005t Massard i et al. 201(1 iTucci et all 
l20ll . The red lines show the sum of the contribution of radio-loud AGN and star-forming galaxies when spinning dust emission is 
included to the SED templates. 



